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ABSTRACT 

We use distant Blue Horizontal Branch stars with Galactocentric distances 16 kpc < r < 48 
kpc as kinematic tracers of the Milky Way dark halo. We model the tracer density as an oblate, 
power-law embedded within a spherical power-law potential. Using a distribution function 
method, we estimate the overall power-law potential and the velocity anisotropy of the halo 
tracers. We measure the slope of the potential to be 7 ^ 0.4 and the overall mass within 50 
kpc is ~ 4 X VS^^Mq. The tracer velocity anisotropy is radially biased with /3 ~ 0.5, which 
is in good agreement with local solar neighbourhood studies. Our results provide an accurate 
outer circular velocity profile for the Milky Way and suggest a relatively high concentration 
dark matter halo (Cvir ^ 20). 
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1 INTRODUCTION 

The mass of our Galaxy is a fundamental - yet poorly constrained 
- astrophysical quantity. Several attempts have been made to mea- 
sure the_Jotal_jTiass_ofttie_Mil^^ kinematic tracers 



Galaxy: kinematics and dynamics — dark matter — Stars 



years has a consensus on th e profile been reached (e.g. lSesar et al.l 
l201lUDeason et al.ll2011bh . Our knowledge of the orbital proper- 
ties of the stellar halo stars is limited to the solar n eighbourhood 



(e.g.lWilkinson & Evansll999l : IXue et alj2008l : lGnedin et al.l2010l : 



where the velocity ellipsoid is radially b iased (e.g. iKeplev et al 
l2007l: ISmith et al.io09iJ - rBond et alj|201oh . In contrast. ISirko et al 



Watkins et al.l |2010|). the orbits of the M agellanic Clouds (e.g. 
Lin & Lvnden-Bel]|l982l : lBesla et alj2"007l) . th e local escape spee d 
dSmith et alj|2007h . and the timing argument jLi & Whitj|2008l) . 
The results of this extensive list of work is distressingly inconclu- 
sive with total masses ranging from 0.5 — 3 x IO^^Mq. 

The most common method to probe the mass distribution is to 
use kinematic tracers such as globular clusters, stellar halo stars 
and satellite galaxies. The properties of these tracer populations 
are linked to the underlying matter distribution via the steady state 
(spherical) Jeans equation: 



dlnptr 
dlnr 



dlncTr 
dlnr 



(1) 



( l2004h inferred the velocity anisotropy of stellar halo stars at larger 
distances (i.e. r > 10 kpc) from line-of-sight velocities alone, 
and found an isotropic velocity ellipsoid; this is in contrast to the 
strongly radial anisotropy found locally. However, the uncertainties 
in these measurements proved too large for a conclusive result. 

In this Letter, we break the mass-anisotropy-density degen- 
eracy for the first time. We adopt the recently measured stellar 
halo density o f iDeason et al] ( l201 Ibl) and disentangle the remaining 
mass-anisotropy degeneracy using line-of-sight velocities of Blue 
Horizontal Branch tracers out to r ~ 50 kpc selected from the 
Sloan Digital Sky Survey (SDSS). Several studies have measured 
the total mass of the Galaxy within r ~ 50 kpc (e.g. iKochanek 
I996I : [Wilkinson & Evanslll999l : ISakamoto et al]|2003l : IXue et al 



At face value, this equation is remarkably simple; the mass distri- 
bution is related to the logarithmic gradients of the radial velocity 
dispersion ct,. and density (ptr) of the tracers, as well as the veloc- 
ity anisotropy (/3). Without firm knowledge of the tracer proper- 
ties, our Galactic mass measures suffer from the well-known mass- 
anisotropy-density degeneracy. 

The density distrib ution of the stell ar halo has been 
studied e xtensively (e.g. Yannv et al. 2000l ; IChen et al] I2OOII ; 



lNewberg & Yannv.2006;iJuric et al .l2008r However, only in recent 
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200^ but the mass profile is poorly known. We measure the slope 
of the overall potential and thus provide an accurate circular veloc- 
ity profile out to 50 kpc. 



2 METHOD 

Beyond the solar neighbourhood, we typically have full spatial in- 
formation for stellar halo stars together with accurate line-of-sight 
velocities. While the density distribution of stellar halo stars has 
been extensively studied, there have been very few attempts to infer 
the velocity anisotropy. In the absence of information on the proper 
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motions, most previous studies have resorted to assuming a velocity 
anisotropy for tiie tracers. Often, this assumption is motivated by 
the predictions of c osmological simulations (e.g. IXue et al ] |2008l : 
iGnedin et allboiol) . However, the line-of-sight velocity distribu- 
tion (LOSVD) itself contains valuable kinematic evidence that is 
often not suitably exploited. For example, information on the ve- 
locity anisotropy is encoded within the fourth order moments of th e 
LOSVD (i.e. the kurtosis; see e.g. Figure 3 of lOeason et al .12011 ah . 
There are two main requirements needed to extract such informa- 
tion from the LOSVD: i) a large sample of tracers and ii) tracers 
with a wide sky coverage. 
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2.1 Halo Tracers: Blue Horizontal Branch Stars 

The Sloan Digital Sky Survey (SDSS) has now mapped an im- 
pressive 20,000 deg^ of sky. Furthermore, the dedicated spectro- 
scopic Sloan Extension for Galactic Understanding and Explo- 
ration (SEGUE) survey has unearthed several thousand kinematic 
tracers of the stellar halo. In particular, ~ 4000 Blue Horizontal 
Branch (BH B) stars have be en spectroscopically identified out to 
r < 60 kpc jXue et alllloT l|). BHB stars are excellent halo tracers 
owing to their intrinsic brightness (Mg ~ 0.5) and accurate dis- 
tance estimates {AMg ~ 0.15). The combination of wide sky cov- 
erage plus a large sample of distant halo tracers provides a unique 
opportunity to break the mass-anisotropy-density degeneracy in the 
halo. In fact, the density p r ofile of these tracers has recently been 
measured bv lDeason et al.l ( 1201 Ibl) (hereafter, DBF) using the lat- 
est photometric SDSS DR8 release. They found that the stellar halo 
out to r ---^ 50 kpc is well-described by a smooth flattened density 
distribution of broken power-law form: 



r-q < 27 kpc, 
r, > 27 kpc, 



(2) 



with Tg — + /q'^, where the minor-axis to major-axis ratio is 
q = 0.59. Armed with this independent measure of the density pro- 
file, a large kinematic sample of BHB stars can be used to constrain 
their velocity anisotropy and the underlying mass distribution. To 
this end, we use BHB stars selecte d from the SDSS DR8 spectro- 
scopic survey bvlXue et al.l 1 I2OI ll) . We assign distances using the 
colour-absolute magnitude relation derived in DBF. For simplicity, 
we only consider stars beyond the ellipsoidal break radius, rq = 27 
kpc (corresponding to spherical radii r > 16 — 26 kpc). This en- 
sures that the tracers are far away from the disc of the Galaxy, and 
allows us to model the density distribution of the tracers with a sin- 
gle power-law. The final sample consists of 1933 stars with Galac- 
tocentric radii in the range 16 < r/kpc < 48. Observed heliocen- 
tric velocities are converted to Galactocentric ones by assuming a 
circular speed0of 240 km s^^ at the position of the sun (Ro = 8.5 
kpc) with a solar peculiar motion {U,V,W)=(ll.l, 12.24, 7.25) 
kms^^. Here, U is directed toward the Galactic centre, V is posi- 
tive in the direction of Galactic rotation and W is positive towards 
the North Galactic Pole. In Fig.[T] we show the final sample of 1933 
BHB stars in the (x, z) plane (left panel) and their line-of-sight ve- 
locities as a function of radius (right panel). 



^ Note that we adopt the recently revised LSR (e.g. iReid et al] |2009| : 
lMcMillaij[20Tlh but our main results are unchanged if we use the conven- 
tional value of 220 kms~^. 



Figure 1. Left panel: The spatial distribution of BHB stais in the x-z plane 
with Tq > 27 kpc. Right panel: The hne-of-sight velocities as a function of 
Galactocentric radius. 
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Table 1. The results of the maximum likelihood analysis. The tracer density 
parameters {a, q) are fixed. 



2.2 Distribution Functions for Splierical and Flattened 
Tracers 

We model the dynamical properties of th e BHB tracers using a dis- 
tribution function (DF) method (see e.g. lBinnev & Tremainell2008l 
and references therein). The phase-space structure of dynamical 
tracers is described by a probability density function; this is often a 
more practical approach than following individual orbits. 

We always assume that the overall potential is spherically 
symmetric, as suggested by several recent studies ( e . g . ISmith et aU 
l009bllKoposov et alJl201(il . lAgnello & Evansll2012h . As we are re- 
stricting attention to stars beyond 16 kpc from the Galactic cen- 
tre, we can safely ignore any flattening influence provided by the 
Galactic disc. For simplicity, we use a power-law profile for the 
potential $ = $0 (f'/lkpc)"''', where 7 is constant. 

We investigate models in which the stellar halo density is a 
spherical power-law, namely p oc r~° with 7 constant. The DFs 
have been described elsewh ere (e.g. lEvans et al .ll997l;lAn&Evansl 
l2006l : lDeasonetalJl2011ii) . The velocity distribution is given in 



terms of the binding energy E 
the total angular momentum L - 



+ + t;D and 



(3) 



(4) 



F{E,L)<xL-'''f{E) 

where 

Here, P is the Binney anisotropy parameter defined as /3 = 1 — 
+ ('^0))/(^r)> and taken as constant. 
We also investigate models in which the stellar halo density 
is a power-law, but with constant fl attening q, as sug gested by a 
number of recent studies (e.g. DBE. ISesar et al.|[201 ll) . The exten- 
sion of eqns ll3}-(|4ll into the flattened regime is given by (see e.g.. 
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Figure 2. The likelihood contour levels for the BHB star's assuming the 
DBE density profile. The blue shaded regions show the Icr (68%) confi- 
dence region, whilst the soHd lines encompass the la (95%) confidence 
region. The gray arrows indicate the confidence regions when a spherical 
tracer density is assumed with slope of 3.5. The assumption of a shallower 
tracer density profile leads to a bias towards lower masses. 



IdeBiuiineet all 19961) 

F{E,L^,lI/L^) (xL-^^f{E)h{e^Ll/L^ 

where 

= 2F1 (1, -; -;e LJL 



(5) 



(6) 



Here, e — \/\ — is the eccentricity, is Pochhammer's 

symbol where (x)k ~ T{x + k)/T[x) and 2F1 is a hypergeomet- 
ric function. Remarkably, the effect of flattening is provided by a 
simple multiplicative factor. 

In our analysis, we assume that the tracer density, described 
by a and q, is known. Our favoured model is the single-power law 
for rq > 27 kpc found by DBE, in which a = 4.6 and q — 0.59. 
However, we also consider the spherical limit of this density (i.e. 
q — l,a — 4.6) and the commonly used spherical density with 
a power-law slope of a = 3.5. These latter two examples illus- 
trate the biases caused when an incorrect tracer density is adopted. 
Note that we ignore rotation in this analysis as previous work has 
infer red that the overall ro tation signal is negligible for the stellar 
halo ( lDeasonetal]|2011al) . 



3 RESULTS 

Our aim is to constrain the overall potential (defined by the slope 
and normalisation) and tracer velocity anisotropy. We have full 3D 
spatial information for the tracers, but only have one velocity com- 
ponent. The LOSVD is constructed by marginalising over the un- 
known tangential components and a maximum-likelihood method 



Figure 3. The line-of-sight velocity dispersion profile of the maximum like- 
lihood models when a DBE density distribution is assumed (blue points) 
and when a spherical power-law profile with index 3.5 is assumed (gray 
points, off set by -50 kms~^ ). The dashed line shows the fit to the BHB data 
derived bv lXue et al.l j2Q08l) . In each case, 2000 stars are drawn at random 
from the maximum hkelihood function. The same hne-of-sight velocity dis- 
persion profile can be reproduced by either model, even though they have 
very difference mass profiles. 



is used to derive the unknown parameters, <l?o> 7 and , 



Here, Fios is the LOSVD (see e.g. equation 7 in iDeason et al.l 
l2011ah and A'^ is the total number of BHB star tracers in our sam- 
ple. The likelihood confidence contours are shown in Fig. |2] when 
our favoured (a = 4.6, q — 0.59) tracer density is adopted. The 
blue shaded region indicates the la confidence region whilst the 
solid black line gives the 2a confidence region. For comparison, 
the gray error bars show the 1 a confidence region when a spher- 
ical tracer density is assumed with slope 3.5, which is a widely 
assum ed value in the literature (e.g. IXue et aflllOOSl : iDeason et al.l 
l2011al) . Our results favour a radially biased velocity anisotropy 
with /3 — 0.5^ ; this is the most accurate measure of the velocity 
anisotropy beyond the solar neighbourhood from line-of-sight ve- 
locities alone. Somewhat surprisingly, our measurement is in good 
agreement with l ocal solar n eighbourhood constraints o n the veloc- 
ity el lipsoid (e.g. lKeplevet al. 2007; Smith et al. 009al: lBond et al.l 
|201(]|) . This suggests that the velocity anisotropy of stellar halo stars 
may be approximately constant over a large radial range. In addi- 
tion, we note that a radially biased velocity anisotropy of /? ~ 0.5 
is in good agreement with the predictions of cosmologic al simu- 
lation s (e.g. Diemand et al. 2007; Sales et al. 2007; Nava rroet al.l 
I2OIOI) . Our results are also co nsistent (within th e errors) with the 
velocity ellipsoid measured bv lSirkoetaLl ( l2004 . While these au- 
thors favour an isotropic ellipsoid, the errors in the tangential com- 
ponents are large and also allow for radial anisotropy. 

This methodology also allows us to measure simultaneously 
the normalisation and slope of the potential, the latter of which has 
not been previously measured for our Galaxy. We measure the over- 
all slope of the potential to be 7 — 0.4 
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Figure 4. The circular velocity profile of the Galaxy. The blue shaded region shows the la constraint found in this work for our favoured tracer density profile 
of DBE. The blue vertical and horizontal line-filled regions indicate the additional uncertainty with systematic errors in the tracer density power-law index of 
±0.2 dex. The gray line-filled region shows the profile when instead a spheiical tracer density with power-law index 3.5 is used. The solid a nd dotted lines 
indicate the maximum likeliho od solutions. Con straints on the circular veloci ty from other studies are sh own by the black eiTor bars:. M cMillan & Binnevl2oTol 
(MB IO; Solar neighbour hood) lMcMillai1 '2011 (Mil: Solar neighbourhood) fwilkinson & Evans|[T99^ (WE99: r = 50 kpc). IXue et al.ll2Q08l (X08: r = 60 
kpc), |GnedinetalJ|2010| (GIO; r = 80 kpc), Watkins et al. 2010 (WIO; r = 100 kpc). The solid/dashed red and green curves show models with dark matter 
components of NFW form and a baryonic component consisting of an exponential disk with mass 5 X 10^" A/q and scale length 3 kpc and a Hemquist bulge 
with mass 5 X 10^ Mq. A dark matter component with virial mass Mvir ~ IO^^Mq and concentration Cvir ~ 20 is favoured by the results of this study. 



within r- = 50 kpc is A/(50) = 4.2±0.4x lO" Mq; in good agree- 
ment with independent mass measures out to this radius (Kochanek 
1 19961 : M(50) = 4.9tl;i x 1 0"Afo[wilkinson & E vans 1999: 
M(50) = 5.4l° j X lO^Mrr,, ISakamoto et aLll2003l : M(50) = 
5.5toi X 10"Mr^.JSmithet alJl2007l : M(5Q) = 3.6 - 4.0 x 
lO^Mrn. lBesla et all2007l : M(50) ~ 4.5 x IO^Mq). 

In Table [T] we give the maximum likelihood solutions when 
different tracer density profiles are adopted. The differences are 
modest when a spherical tracer population is assumed; the velocity 
anisotropy is less radially biased and the potential slope is slightly 
shallower. On the other hand, adopting a shallower tracer density 
slope (i.e. a = 3.5) leads to significant differences in the total 
mass; the mass is 30% higher when the tracer density profile of our 
favoured model (a — 4.6) is assumed instead of the more com- 
monly adopted a — 3.5. 

In Fig. (3) we show the line-of-sight velocity dispersion profile 
of our maximum likelihood solutions. We draw 2000 stars from the 
appropriate distribution function using a Monte-Carlo method; the 
blue and gray points show the resulting line-of-sight velocity dis- 
persions for the DBE (a = 4.6, q = 0.59) and spherical compari- 
son (a = 3.5) density models respectively. The gray points are off- 
set by -50 kms^^ for illustrative purposes. The dashed l ine shows 
a fit to the observed relation for BHB stars derived by IXue et all 
( l2008h . This plot highlights the degeneracy between different den- 
sity models; both are able to reproduce the observed line-of-sight 
velocity distribution but they give very different mass profiles. 



3.1 Milky Way Circular Velocity Profile 

In Fig. [4] we show the circular velocity profile of the Galaxy de- 
rived from our model. The blue shaded region shows the la confi- 
dence region from this work assuming a DBE density profile with 
slope a — 4.6 and flattening q — 0.5^ We also consider sys- 
tematic uncertainties of 0.2 dex in the tracer density slope. Note 
that we do not consider systematic uncertainties in the flattening of 
the tracer density profile, as the mass profile is much more sensi- 
tive to changes in the power-law index of the tracer density (e.g. 
see Table [TJ. DBE found statistical errors of ~ 0.1 — 0.2 dex in 
the power-law indices and systematic effects due to the presence 
of un-relaxed substructure can also cause biases of ^ 0.2 — 0.3 
dex. The blue line-filled regions indicate the additional uncertain- 
ties that could be caused by such biases. The normalisation of the 
circular velocity is slightly increased or decreased when the tracer 
density is modified upwards or downwards by 0.2 dex. However, 
the overall potential slope and tracer velocity anisotropy are hardly 
changed. This emphasises that an accurate measure of the tracer 
density profile is vital to infer the mass profile of our Galaxy. The 
hatched gray region gives the mass profile when a spherical den- 
sity profile with slope a — 3.5 is assumed instead. The points with 
error bars show other constraints from the literature. 

The green and red lines show model profiles assuming a bary- 
onic component consisting of a (spherically averaged) exponen- 

^ We note that our assumption of a single power-law for the overall poten- 
tial is most uncertain at the end points of our radial range (i.e. r < 20 kpc 
or r > 45 kpc). 



© 2012 RAS, MNRAS 000,[T]l5] 



Broken degeneracies 5 



tial disc with mass 5 x IQ^^Mq and s cale length 3 kpc and 
Hemquist bulge wi th mass 5 x 10^ Mr (cf.lKlvpin et alj2002l and 



iGnedin et alj|20ld : see also lSinnev & Tremai 



2008) and a dark 



matter component of Navarro-Frenk- White ( Navarro et al.|[l996l) 
fom. The red and green lines assume virial masses of Mvir = 
1 X 10^^ Mq and Afvir = 2 x 10^^ M© respectively, while the 
solid and dashed lines are for halo concentrations of Cvir = 20 and 
Cvir = 10. 

An NFW model with a virial mass of A/vir ~ IO^^A/q and a 
relatively high concentration, Cvir = 20 (solid, red line) is favoured 
by our new constraint on the Galactic potential in the radial range 
16 < r/kpc < 48. A virial mass of 10^^ Mq is towards the 
lower end of independent me asurements from local group timing 
arguments jLi & White'2008h and the kinematics of satellite galax- 
ies lWatkins et aU (i2010i). However, constraints from kinem atic stel- 
lar halo tracers (e.g.lBattaglia et al .120051: IXue et al.l2008h . the local 
escape speed /Smith et al| 2007 ) and the orbits of the Magellanic 
clouds (e.g Busha et al. 2 01 ih also favour a less massive halo with 
A'/vir ~ 0.5 - 2 X IO^^A'/q. 

The mass-concentr ation relation of cosmological simulations 
(e.g iMaccio et alj|2008t) predict a mean concentration of Cvir ~ 
10 for haloes of mass 10^^ A/0. Thus, the concentration of our 
favoured dark matter halo seems to be at odds with the predictions 
of ACDM simulations. However, the inferred mass-concentration 
relat ion is based on dark ma t ter only simulat i ons. Several authors 
(e.g iBlumenthal et alj [19861 : IMo et all 1 19981 : iGnedin et al.ll2004h 
have suggested that the influence of baryons at the centre of the 
halo potential well can lead to an adiabatically contracted dark mat- 
ter halo (i.e. a more high l y conc entrated dark matter component). 
We note that I Smith et alj ( l2007h favour a adiabatically contracted 
dark matter halo from their local escape velocity constraints. The 
'standard' dark matter halo of this work is strikingly similar to the 
solid, red model shown in Fig. |4] with Afvir = 0.9 x lO'^^ and 
Cvir = 24. Note that the concentration of the favoured model of 
ISattaglia et alj ^20051) is also relatively high for their inferod halo 



A/v 



■ 0.8 X lO"-" and Cvir 



18. 



4 CONCLUSIONS 

We studied the potential of the Milky Way by using BHB stars as 
halo tracers. We modelled the tracer density as a flattened power- 
law embedded in a spherically symmetric, power-law gravitational 
potential. The shape of the tracer velocity distribution is controlled 
by the constant velocity anisotropy parameter /3. We used a maxi- 
mum likelihood method to derive the potential power-law slope 7 
and normalisation $0, in addition to the velocity anisotropy of the 
tracer population. 

By adopting our favoured tracer density model with slope 
Q = 4.6 and flattening q — 0.59, we find that BHB tracers in 
the radial range 16 < r/kpc < 48 have a radially biased velocity 
anisotropy with /? = 0.51q'2**. The agreement of this result with 
local solar neighbourhood measurements suggests that the veloc- 
ity anisotropy of stellar halo stars may be constant to a good ap- 
proximation over a relatively large radial range. We also measure 
the power-law slope of the overall potential to be 7 ~ 0.4, which 
lies in-between the isothermal (7 — 0) and Keplerian (7 = 1) 
regimes. This model implies that the total mass within 50 kpc is 



^ Note that smaller scale lengths (e.g. 2 kpc) make little difference to the 
circular velocity profile in the radial range probed by this work. 



A/(50) ~ 4 X IO^^A/q, in good agreement with other independent 
estimates in the literature. Neglecting the flattening of the stellar 
halo only has a small effect; the velocity anisotropy is less radially 
biased and the slope of the potential is slightly shallower. However, 
shallower tracer density profiles - such as the commonly adopted 
power-law a — 3.5 - lead to a lower inferred halo mass and a 
different circular velocity profile. 

For the first time, we have provided a measure of the circular 
velocity profile of the Galaxy in the Galactocentric radial range 
16 kpc < r < 48 kpc. Our results suggest that the dark matter 
potential may be more centrally concentrated than the predictions 
of dark matter only simulations; this might be a consequence of 
adiabatic contraction in this inner radial regime. 

This work is a useful step towards obtaining a detailed descrip- 
tion of the mass profile of our Galaxy. However, further progress 
requires tighter constraints on the tracer population properties at 
larger distances (r > 80 kpc). This presents a daunting task; not 
only do we suffer from a lack of kinematic tracers st such large dis- 
tances, but we have very little knowledge of the tracer density pro- 
file. Moreover, we have no constraint on the velocity anisotropy as 
we purely observe the radial velocity component at such distances 
(i.e. vios ~ Vr). New and future surveys and observational facili- 
ties, such as the Large Synoptic Survey Telescope and the planned 
30m class of telescopes, will be of vital importance in order to 
tackle this problem. 
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